Existence results for integral viscoelastic fluids 

Laurent Chupin0 
Abstract 

We consider the flows of viscoelastic fluid which obey a constitutive law of integral type. The existence 
and uniqueness results for solutions of the initial boundary value problem are proved, and the stationary 
case is studied. 



1 Introduction 

The objective of this paper is to provide mathematical results on the integral models for viscoelastic flows: 
existence solution, uniqueness and stationary study. 
• ■ In integral models, stress components r are obtained by integrating appropriate functions, representing the 
I , amount of deformation, over the strain history of the fluid. The constitutive law given the stress at time t and 
at position x is written as follows: 

<N . t 

>; T{t,x)^ m{t-T)S{F{T,t,x))dT. (1) 

t The scalar function m (the memory) and the tensorial function S are given by the properties of the fluids 
' studied, whereas the deformation tensor F is coupled with the velocity field of the flow. This flow is itself 
, governed by the Navier-Stokes equations, this constitutes a very strong coupling between the velocity and the 
' stress. 

p ^ , In some rare cases, it is possible to express the integral models into differential forms (as in the Maxwell models 
^rfH which corresponds to the case where S is linear and where the memory m exponentially decreases). For these 
. , differential models, there are many mathematical results in the same spirit of those presented here (see for 
^ ■ instance UHl [HI HI [H [Ml [3D] ) ■ 

, But for really integral models there are far fewer results. One of the only relevant work on this type of model 
^ ' is that of M. Renardy [33]. In its paper, M. Renardy prove an existence and uniqueness result for a K-BKZ 
I— I, fluid using Kato's theory of quasilinear hyperbolic equations. Its elegant approach differs substantially from 
the approach used here, and does not seem easily adaptable to more general laws. For instance, in [33] the 
, author assumes that m is not singular at contrary to what is predicted by some molecular models like the 

■ Doi-Edwards model. It is important to note that while the theoretical results are very few, many authors have 
P\| , studied the numerical simulation of flow with an integral law of type (dJ . The review article [2S] and references 

■ cited therein, provide a good overview of the state of the art regarding the various methods. 

, The main reason for this lack of theoretical results is probably due to the nature of the equations: 
T-H ■ • To evaluate the stress in an integral model, we must know all the previous configurations. This difficulty 
, is overcome by introducing an additional time variable corresponding to the age. It is then necessary to 

■ manage two different times, and in particular use a Gronwall type lemma in two variables in order to obtain 
. . , fine estimates with respect to these variables. 

^ ' • The usual integral models are strongly nonlinear (in the linear case we find the well-known Maxwell 
k> , models). The possibility to circumvent this difficulty is to work with solutions regular enough, more precisely 
^ ■ in a Banach algebra like the Sobolev spaces W^'P for p large enough, typically p greater than the dimension 
of the physical fluid domain. 

Organization of the paper — Section [2| is devoted to the presentation of the model. The dimensionless 
form of equations, as well as many classic examples are given. The main results are stated in the Section [3] 
whereas the proofs are given in the next sections. Section[3|is entirely devoted to the proof of the first theorem 
regarding the local (in time) well posseness. The three next sections ([SJ \E\ and [7]) concern the proof of the 
uniqueness result, the global existence with small data, and the case of stationary solutions respectively. The 
conclusion of this paper (section [8]) contains many remarks and open questions. Finally, some notions on 
tensors, and a technical Gronwall type lemma has been postponed to Appendices [XI and [B] 
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2 Governing equations 



2.1 Conservation principles 

The fluid flows is modeled using the equation of conservation of the linear momentum and the equation of the 
conservation of mass, which read in the incompressible and isothermal case as follows: 

p {dfU + u ■ Vm) + Vp = div cr + f, 
divM = 0, 

The real p is a constant density of mass and the vector / corresponds to some external body forces. This 
system is closed using a constitutive equation connecting the stress and the deformation Du — i(VM+ '^{Wu)). 
For a so-called Newtonian viscous fluid, the relationship is linear: Tg — 2rjsDu. The real 77^ > is named the 
solvent viscosity and the contribution divx^ in the momentum equation gives the usual diffusive term 77^ Am. 
To taking account some elasticity aspect, appearing for instance in polymer solution, we add to the viscous 
contribution an elastic one: 

(T = 2-qsDu + r. (3) 

The role of this additional contribution t is to take into account the past history of the fluid. The most natural 
way to do this is to introduce the integral models. 



2.2 Integral models 

Very generally, the elastic contribution r(t, x) at time t and at spatial position x is written 

T{t,x)= S {F{T,t,x)), (4) 

where 5^ is a functional to clarify, which depends on the deformation gradient F{T,t, •) from a times T to a 
next time t. 

More precisely, the deformation gradient F{T,t,-) measures stretch and rotation. It is deflned as follows: 
for two times T < t given, we first introduce the notation x(T, t, X) which corresponds to the position at 
time t of the fluid particle which was at the position X at time T. The dynamics of any mechanical problem 
with a velocity field u(t, x) can be described by this flow map x{T, t, X) which is a time dependent family of 
orientation preserving diffeomorphisms: 

idtx{T,t,X) = u{t,xiT,t,X)), 

I x{T,T,X) =X. ^ ' 

The deformation gradient F(T, t, X) is used to describe the changing of any configuration, amplification or 
pattern during the dynamical process. It is defined by 

FiT,t,X) = —{T,t,X). (6) 
Finally, the deformation gradient F{T,t,x) will be defined as the corresponding in the Eulerian coordinates: 

FiT,t,x{T,t,X))^F(T,t,X). (7) 
The integral models we study in this article correspond to the particular case of equation They are written 



T{t,x)^f m{t-T)S{F{T,t,x))dT, 

J —00 



(8) 



where m is called a memory function and 5 is a model-dependent strain measure. It is a tensorial function 
(its arguments and its images are 2-tensors). 

Remark 2.1 Due to many physical principle, the functions m and S must satisfy some assumptions. 

• For instance, the principle of frame indifference implies that the stress tensor depends on the relative defor- 
mation gradient F only through the relative Finger tensor '^F ■ F (or its inverse, the Cauchy-Green tensor), 
see the examples given in the subsection \2.4\ 

• In the same way (see also the subsection \2.4\ l, the principle of fading memory implies that m must be a 
positive function which decreases to 0. 
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2.3 A closed system 

2.3.1 PDE for the deformation gradient 

We derivate the relation ^ with respect to the time t. The chain rule together with the relation ([5]) yields 
the following equation 

dtF{T, t, X) = dtF{T, t, x) + dtx{T, t, X) ■ d^F{T, t, x) 
^ dtF{T,t,x) +u{t,x) ■ d^F{T,t,x). 

But using the relation © together with the chain rule and the relation ([5]) again, we obtain 

dtF{T,t,X)=dx{dtx{T,t,X)) 

= dx{uit,x{T,t,X))) 
^dxx{T,t,X)-dMt,x) 
= F{T,t,x) ■ d^u{t,x). 

The equations ([9]) and ((TO)) show that we have the following relation coupling the velocity field u and the 
deformation gradient F: 

dtF + U- VF = F -Vu. (11) 

2.3.2 A new time variable to take into account the past 

Note that in the previous subsection, the time T can be view as a parameter. In fact, it is only used in the 
law Q, or in the law ([5]) for the integral form, as a marker of past events. In the sequel, it is interesting to 
select as independent variable the age s = t — T, which is measured relative to the current time t. We now 
introduce G{s, t, x) = F{t — s, t, x). Clearly, we have the following relation instead of the relation ([TT|): 

dtG + dsG + u-WG = G-Wu, (12) 

where naturally the velocity u only depends on (t, x), and is independent of this new variable s. Moreover, in 
term of variables {s,t), the relation (|8]) given the stress tensor reads 

r+co 

T{t,x)^ m{s)S{G{s,t,x))ds. (13) 

Jo 

Initial and past conditions — The equation (TT^ requires initial conditions, that is to say that we must give 
the value G{s,0,x) at time i = for all age s > and for all point x E ft, and the value G{0,t,x) at any 
time t >0 and for all point x E ft. Note that for T fixed, the deformation field F{T,t,x) can be thought of 
as having been created at time t = T with the natural initial condition F{T,T,x) = d. In term of the new 
variables {s,t) this relation reads 

G{0,t,x) = d. (14) 

The condition about G(s, 0, x) correspond to the description of the deformation gradient before the time t = 0. 
In practice this value is unknown but it is physically reasonable to assume that for an age old enough, the 
fiuid was quiescent. In practice we will only assume that there exists a given function Goid such that 

Gis,0,x) ^Gomis,x). (15) 

2.3.3 Non-dimensional final model 

The resulting system using the equations ©, (O, ([T^ and ([T^ is then written 

p {dtu + u ■ Vm) + Vp — Am = div r + f, 



div M = 0, 



T(t,x) — / m(s) 5(G(s, t, a;)) ds, 
Jo 

dtG + dsG + u-VG^G- Vu. 



(16) 
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This system can be adimensionalized in the usual way. We introduce the characteristic values U and L for 
the velocity and the length. The current time is then of order of L/U and it is natural to introduce another 
characteristic time A for the age variable s. The characteristic viscosity of the fluid takes into account the 
viscosity rjs of the solvent, but also the viscosity rje of the elastic part (the polymer): we note rj ^ rjs + rje- 
More precisely, we introduce the following dimensionless variable, quoted by a star: 



t* 



t 



L/U' 



P 



P 



rjU/L' 



r 



f 



T]U/L^ 



" -, G* = G, 5^(G*) = ^, m^is^)^"^^'^ 



(17) 



The system (|16|) is then written in dimensionless form as follows (where we drop the star for sake of simplicity) : 



' 5He(9tM + u ■ Vtt) + Vp - (1 - uj)Au = divr + /, 
div u — 0, 

uj /■+°° 

T(t,a;) = ^y m{s)S{G{s,t,x)) ds, 
dtG + r^dsG + uVG = GVu, 
where we introduced the three non-dimensional numbers which characterize the flow: 



(18) 



The Reynolds number 9le 
acting on the fluid; 



which corresponds to the ratio between inertial and viscous forces 



• The Weissenberg number = which is the ratio between the time of the relaxation of the fluid 
and the time of the experiment; 

• The retardation parameter w = ^ G [0, 1] which balances the purely viscous effects (w — 0) and the 
purely elastic effects [lo = 1). 

The system (fT5| is closed with the following initial and boundary conditions: 



it=o 



Wo, 



0, 



G\ 



Go\i 



G 



s=0 



(19) 



Remark 2.2 As we precised in the Remark \2.1[ the stress tensor t depends on the deformation tensor G via 
the right relative Finger tensor B = ^G ■ G or via the Green-Cauchy tensor C — . Using the last equation 
of (jlSp . we note that the tensor B satisfies 

dtB + -^dsB + uVB ^ B ■Vu+ ^(Vm) • B, (20) 

whereas the tensor C satisfies 

dtC + -^dsC + u VC = -C ■ '^(Vm) - Vm • C. (21) 

2.4 Examples of integral models 

In this section, we present some classical integral laws of kind (IT51) to model the viscoelasticity. These law are 
defined by the memory function m and by the strain measure S. 
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2.4.1 Memory function m 

Usually, in the viscoelastic formulas, a relaxation function g is introduced to describe the return into equilib- 
rium. It physically corresponds to the response of the stress to a shear jump. In accordance with thermody- 
namics through what is called the principle of fading memory - see |11) - this relaxation function is constrained 
in that g(0) = 1 and lim+oo 3 = 0. To so called memory function m used in the integral models (fT3)) corresponds 
to m = —g'. Consequently, the function m is decreasing, positive and satisfies m[s) As = 1. 
In many cases experimentally observed relaxation functions exhibit a stretched exponential decay e~'--'^l^^ where 
A > is a relaxation time. We could as well have considered the case of several relaxation times, that is a 
memory function like 

N 

m{s)=Y.^^'''^''- (22) 

From a mathematical point of view, it will be equivalent to consider m(s) = e^^ (on the dimensionless form). 
This expression for the memory function can be generalized. For instance, in the Doi-Edwards model - see |13) . 
the memory function is given by 

^(s) = -A-^e-(2'=+i)^/\ (23) 

fc=0 



In practice, the difference between this model (j23l) and the model (j22p containing a finite number of relaxation 
times is really important. In fact, in the case of the model (j23p . the function m has a singularity in 0. This 
singularity can bring additional difficulties (eg, such a case is not treated in the article [33] of M. Renardy). 
The function m remains integrable, which is the key assumption for the present results. 

Even if the exponential case is usually used, many other possible choices for the memory function m are 
possibles - see [IS]. The algebraic pattern g{s) = {s/\)~^ with < /3 < 1 is observed in the stress relaxation 
of viscoelastic materials such as critical gels [5] I35j , in the charge carrier transport in amorphous semicon- 
ductors [33] , in dielectric relaxation [JB] or in the attenuation of seismic waves [37] . That corresponds to the 
following memory functions 



k=l 

2.4.2 Strain measure S 

MsLxwell models — The more simple case corresponds to the choice S{G) ^ B 5 where B ~ ■ G, and 
where the memory (dimensionless) function m is given by rn(s) — e^'*. The stress tensor t is then given by 

T{t,x) = —j^ e-^{B{s,t,x)~d)ds. (25) 

This expression is simple enough to deduce a PDE for the stress tensor r from the PDE for the deformation 
tensor G. In fact we use the equation ([30]) satisfied by the Finger tensor: 

dt{B -S) + —ds{B -S)+u- V(B - S) = {B - S) ■ \/u + "^(Vm) ■ {B - 6) + 2Du. (26) 

We next multiply this equation by ^'7i(s) and integrate for s e (0, -l-oo). Taking into account the initial 
condition B\s=o — we obtain the Upper Convected Maxwell (UCM) model: 



OTe 



(^dtT + M • Vr - ^(Vit) - T-T ■ Vv}j + T = 2luDu. (27) 

Another classical case is to the Lower Convected Maxwell (LCM) model. It corresponds to S{G) = 5 — C 
where C = B~^, and to an exponential memory function m(s) = e"". Using the equation (|21[) we obtain, like 
to get the UCM model: 

me(dtT + U-VT + T- ^(Vm) + Vm • + t = 2ujDu. (28) 
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Remark 2.3 In fact there exists a continuum of such model (called Oldroyd models, corresponding to a balance 
between upper- convected and lower- convected) but we do not know if these models derive from integral models. 



K-BKZ models — Among the most relevant non-linear cases, the most popular integral models for a vis- 
coelastic flow are the K-BKZ models introduce by B. Bernstein, E. A. Kearsley and L. J. Zapas jiJIS] and A. 
Kaye [21]. For such models, S takes the following form: 

S{G) = Mh,l2KB-d) + Mh,l2){S-C), (29) 

where 0i and 02 are two scalar functions of the strain invariants Ii = Ti{B) and I2 = Tr{B^^) (see the 
Appendix |2] for a discussion on these invariants) . Clearly, the two Maxwell models presented in the previous 
paragraph are particular K-BKZ models, for {<f)i,<p2) — (1,0) and ((/)i, (/12) — (0, 1). 

Remark 2.4 Add a diagonal tensor to S{G), that is consider S{G) + (f>S instead of S{G), do not change the 
mathematical structure of the system. In fact t becomes t -f and the additional contribution can be 

consider as a pressure contribution. 

Following the last remark, we can see the PSM models presented by A. C. Papanastasiou, L. Scriven and 
C. Macosko in [32 as K-BKZ models: 

S{G) = h{hj2)B with /j(/i,/2) = — — -— -. (30) 

a + ph + (1 - p)h - 3 

In these models, the parameters a > and < /3 < 1 are obtained from the rheological fluid properties. In 
the same way, Wagner [37l [38] proposes the following law: 



S{G) = h{hj2)B with h(hj2) = expi-a^ph + (1 - /3)/2 - 3). (31) 

Doi-Edwards model — The Doi-Edwards model is a molecular model where the motion of the polymers is 
described by reptation in a tube, more precisely it corresponds to the simplest tube model of entangled linear 
polymers. The memory function m associated to such model is given by the relation (j23p whereas the strain 
measure is obtained as an average with respect to the orientation of tube segments. Works of P.-K. Currie 
show that we can approach this model using the following strain function (named the Currie approximation, 
see [12): 

S(G) = — ^ — -B - — ^ ^, C with J = h+ 2JI2+3.25. (32) 

^ ' 3(J-1) 3(J- l)V/2 +3.25 1 V 2 y J 



3 Main results 

3.1 Mathematical framework 

In the sequel, the fluid domain $7 C K*^, d > 2, is a bounded connected open set with a Lipschitz continuous 
boundary d^. We use the following standard notations: 

• For all real s > and all integer p > I, the set W''''P{rt) corresponds to the usual Sobolev spaces. We 
classically note LP{n) = W°^P{n) and iJ^(17) = W^^in). 

We will frequently use functions with values in or in the space /^(R*^) of real d x d matrices. In all cases, 
the notations will be abbreviated. For instance, the space {W^'P{il.))^ will be denoted W^'P{il.). Moreover, all 
the norms will be denoted by index, for instance like wi,p(n)- 

• Since we are interested in the incompressible flows, we introduce IIp{i}) ^ {v e LP{V,) ; divv — 0, v n — 
on dfl}, where n is the unitary vector normal to dU, oriented towards the exterior of ft. Moreover, we 
note V{n) = H2{n) n H^{n) and V'{n) its dual. 

• The Stokes operator Ap in Hp{fl) is introduced, with domain D{Ap){fl) = W'^'P{fl) n W^'^{fl) Ci IIp{n), 
whereas we note D^p{n) ^ {v e Hp - ||v||Lp(n) + ( /o^°° WApQ-^^^vW^^,^. dtf'' < +00}. 
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• The notation of kind L^{0,T; D{Ap)) denotes the space of r-integrable functions on (0,T), T > 0, with 
values in D{Ap). Similarly, expressions like g G L°°(R+; L''(0, T; _L''(f2))) means that 

rT 1 

sup ( / \\g{s,t,-)\\lpr^.dt)' < +00. 

Assumptions — About the model ([T5| itself, it uses the two given functions m and S. From a mathematical 
point of view we assume very general assumptions (satisfied by all the physical models introduced earlier) : 

^+00 

(Al) TO : s e R+ I — > m{s) e M is measurable, positive and satisfies / to(s) ds — 1; 

Jo 

(A2) S -.Ge CiR"^) ^ S{G) e CiR"^) is of class C\ 

Note that the notion of derivative for the 2-tensorial application S will be precised in the Appendix El 
Moreover, if we wanted to be more precise, the assumption (A2) is written rather "the function S is of class 
on a subset of £(R'') taking into account the fact that detG = 1" - see Appendix 1X1 again. Throughout the 
remainder of this paper, these two hypotheses (Al) and (A2) will be assumed satisfied. 

3.2 Statements of main results 

The first result concerns an existence result for strong solutions. It is obviously local with respect to time (as 
for the results on the Navier-Stokes equations): 

Theorem 3.1 (local existence) Let T > 0, r +oo[ and p E]d, +oo[. 

Ifuo e Dl,{n), Gold e L°°iR+;W^^P{n)), a^Goid e i''(K+;£P(l^)) and f G L'^(0, T; LP(rj)) then there exists 
Ti, £]0,r] and a strong solution {u,p,T,G) to the system (jlSp in [0,Ti,], which satisfies the initial/boundary 
conditions (jl9p . Moreover we have 

ue U-iO,n;W^'P{n)), dtueUiO,T,;LP{n)), 

TeL°°(0,T,;VKi'f(r!)), 9tT e L'-(0,r,;i?'(r!)), (33) 

G e L°°(K+ X (0, T,); W^^Pifl)), d,G, dtG e L°^{m+;L^{0, T,; LP{n))). 

We will show that the solution obtained in the theorem 13.11 is the only one in the class of regular solution. 
Precisely, the result reads as follow 

Theorem 3.2 (uniqueness) Let T > 0. 

Ifua e H, Gold e L°°{R+;L^{n)), f e L^{0,T;V' {ft)) and if possess two weak solutions (mi,pi, n, Gi) 

and (m2,P2, ''"2, G2) in the usual sense, with for i g {1,2}, 

u, e i°°(0, T; L^{n)) n L^{0, T; H^{il)) D L^{0, T; W^'°°{n)) 

G, e L°°(R+x(0,T);H/i^'^(17))), ^^^^ 
then they coincide (pi and p2 coincide up to an additive function only depending on t). 

If the data are small, it is possible to show that the unique local solution obtained by the theorem 13.11 and 1 3.21 
is defined for all time (up to an assumption on the relaxation parameter a;): 

Theorem 3.3 (global existence with small data) Let r e]l, +oo[ and p ejd, +oo[. 

For each T > 0, with the same assumptions as in the theorem \3.1l there exists ujt G (0, 1) such that if 
< oj < LOT o,nd if the data uq and f have sufficiently small norms in their respective spaces, then there exists 
a unique strong solution {u,p,T,G) to system ([T5)) - (fro)) in [0,T] which belongs in the same spaces that the 
local solution obtained in the theorem \3.1\ 

The last theorem which is proved in this paper concerns the stationary solutions, that is solutions which do 
not depend on time. Since the mathematical problem is fundamentally different, the proof that we propose 
(see Section [7]) requires some additional assumptions on the memory function m and on the function S. 
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(A3) There exists a > and c > such that m(s) < ce for all s G IR+. 

(A4) The function S and its derivative have polynomial growth: there exists (a, b, c) G M'^ such that for all 
G e we have \S{G)\ < c|G|'^ and \S'{G)\ < c\G\\ 

Note that the assumption (A3) implies that the memory m decays sufficiently fast (exponentially) at infinity, 
and prohibits a singularity at 0. The assumption (A4) is generally satisfied by classical law (see the Appendix 1X1 
in which we show that the assumption (A4) is satisfied with a = and 6 = —1 for a PSM law). 
Under these additional assumptions (A3) and (A4), we have: 

Theorem 3.4 (stationary solutions) Let p G]d, +oo[. 

If Lj is small and ?// G L''(f2) have a sufficiently small norm then there exists exactly one small strong solution 
{u,p^T,G) to the stationary problem associated to (jl8p - (|19p with 

ueW^-Pin), reW^^Pin), G e L°°{R+{d^i);W^^Pin)), dsG e L°°iR+{d^i)-LP{n)), (35) 

where d/i ~ e^^'^ds, (3 being a positive number (depending on p, d, uj and Wt). 

Remark 3.1 Here are some remarks concerning the four above theorems: 

- We remark that the theorem \3.3\ do not contain smallness condition on the deformation tensor Gold- -^t rest 
this tensor is not zero but is equal to the identity tensor. An assumption of smallness should eventually be 
introduced on the quantity Gom ^ ^- In fact it is implicit in the smallness assumption on the parameter uj. In 
the same way, the term "small solution" used in the theorem \3.4\ does not concern the tensor G. 

- The theorem \3.S\ can be viewed as a result of stability for the solution u ~ 0, t ^ 0. 

- The assumptions (A3) and (A4) are not optimal (a more precise formulation will be complex). For instance, 
following the proof of the theorem \3.4\ (see section^ we can see that ifa — b+l = then the assumption (A3) 
is unnecessary. 

- An example of value for a and b is given in the Appendix \^ for the PSM model (more precisely for the 
model (I30p with a = 4 and /3 = 1 j. In this case we have a = & + f = 0. 



4 Proof for the local existence 

This section is devoted to the proof of the local existence theorem l3.1l The main ideas to prove the theorem l3.1 
are based on works of J.C. Saut and C. Guillope [301 HI]. Roughly speaking, we rewrite the equations (1T51) as a 
fixed point equation and applying the Shauder's theorem. This principle was taken up by E. Fernandez-Cara, 
F. Guillen and R.R. Ortega [TH[TS] in the context of the functional spaces U — Lp. That is this choice which 
is presented in the present paper. 

We then first analyze three independent problems. A linear Stokes system with given source term and initial 
value 

<Kc dtu + Vp - (1 - w) Am = g, 

divM = 0, (36) 
u\dn = Q, M|t=o = Mo; 

A following tensor problem given the deformation gradient as a function of a given velocity field m: 

<9fG'+-^a,G + M-VG = G- Vm, , , 

!2ne ' (37) 

G|s=o = ^1 G|i=o ~ Gold; 



And the constitutive integral law given the stress tensor t with respect to a deformation gradient G: 



^ r+<x> 

■(i,a;) = ^y 'm{s)S{G{s,t,x))ds. 



(38) 
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4.1 Estimates for the velocity u solution of a Stokes problem 

The results for the Stokes systera ([5^ are very usual. In this subsection we only recall, without proof (we can 
found a proof in |18|). a well known results for the time dependent Stokes problem: 

Lemma 4.1 Let T > 0, r e]l, +oo[ and p e]l, +oo[. 

If Uq G Dp{fl) and g £ L^(0,T;iJp) then there exists a unique solution u € L^{0,T; D{Ap)) such that 
dtU e {0,T; Hp) to equations (j36|) . Moreover this solution satisfies 

\\u\\L^{f),T-w^^p(Q.)) + \\dtu\\L^{o,T-Lp(n)) < (^e ||Mo||w2.p(n) + ||9||l'-(o, T;Lp(a))), (39) 

where the constant C'l only depends on Q, r and p. 

4.2 Estimates for the deformation gradient G 

This subsection is devoted to the proof of the following lemma, which gives estimates for the solution to the 
system ((37| : 

Lemma 4.2 Let T > 0, r e]l, +oo[ and p e]d, +oo[. 

If Gold e L°^{R+;W^^P{n)), dsGoid e L'' {R+ ; LP (n)) and Tie L''{0,T; D{Ap)) then the problem 1^ admits 
a unique solution G G L°°(R+ x (0, T); M^i'P(r2)) such that dsG,dtG G L°°(R+; i'~(0, T; iP(r2))). Moreover, 
this solution satisfies 

||G||Loo(R+x(o,T);Wi.P(n)) + ||5sG||ioc(R+.2,r(o ^.iP(f2))) + ||9tG||L°=(R+;L'-(0,T;LP(O))) ^^^^ 

< C2(l + ||VM||Lr(o,T;Lp(n)))exp(C3||VM||ii(o,T;vyi'!'(a))), 

where the constants C'2 and C3 depend on f2, p, r, and the norms of Gold md dgGoid- Their expressions 
will be precised in the proof. 

The existence of a unique solution to ([57]) follows from the application of the method of characteristics. Some 
details are given in [15| (Appendix p. 26) about the equation ([TT|). that is using the function F. In practice, 
the following estimates will be made on regular solution G„ which approaches the solution G when a regular 
velocity field m„ approaches the velocity u. The regularity of these solutions G„ with respect to t and s comes 
from the Cauchy-Lipschitz theorem. For sake of simplicity, we omit the indexes "n". In the following proof, 
we refer to |15j for the passage to the limit n — >■ +00. The rest of the proof of the lemma is split into three 
parts: in the first one (see the subsection 14. 2 . 1( 1 we obtain a first estimate concerning the regularity of G, and 
in the subsection 14.2.31 we obtain the estimate for dtG. This estimate requires an estimate on dsG, which is 
given in the subsection 14.2.21 Note that we strongly use a Gronwall type lemma whose the proof is given in 
the Appendix [Bl 

4.2.1 Estimate for the deformation gradient G 

Let p > d and take the scalar product of the equation ((37l) by |G|p~^G. We deduce 

+ + V(|Gn = |Gr^(G . V«) : G. (41) 

Integrating for x £ fl, due to the incompressible condition divU = and the homogeneous boundary Dirichlet 
condition for the velocity, we obtain 

dtm\Un)) + ^^M\\G\\Un))=pJjG^^^^ (42) 

We next use the continuous injection W^'P{n) ^ L°°{n), holds for p > d and making appear a constant 
Co ^ Coin, p): 

dti\\G\\l,,n)) + ^M\\G\\Ui^^^ (43) 
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Now, we take the spatial gradient in ([57]) and compute the scalar product of both sides of the resulting equation 
with jVGI^'^^VG (we will note that this is a scalar product on the 3-tensor, defined hy A :: B = Ai,j^kBij^k)- 
After integrating for x G fl we obtain 

dtiW^GWl^in)) + < iVCnVUl +pj^ \G\\WGr'\W'u\. (44) 

Using the Holder inequality, we have 

dt{\\^G\\l,^,,^) + ^^d.{\\WG\\l,^^^)<2p\^^^^^ (45) 
For p > d, using the continuous injection W^'P{fl) ^ L°°(n) again, we deduce 

5*(l|VG|IL(a)) + < 3pC7o||Vll|U..(^)||Gr,^,.,(^). (46) 

Adding this estimate (|46l) with the estimate (|43p . we deduce 

5t(||G|ki.P(n)) + ^9,(||G||„.i,p(o)) < 3Co||Vm||m/i,p(o)||G||h/i,p(o). (47) 

Using the initial conditions we have ||G||vi/i.p(i2) |^^q = \/d|r2|p and |lG||vyi'P(a) I^^q = llGoidllTyi.p(n), the 
Gronwall type lemma (see the Appendix |B]) implies that for all (s, t) e x (0, T) we have 

||G(s,t,-)||iyi.p(n) <C(s,t)exp(3Co^ II Vm|1m/i,p(o)) , (48) 

whereC(M)4"''°'^""^'^^r'^'^^^ ^f^^^^' 

The assumption Gow e L°°(M+; VKi^P(fi)) implies C e L°°(R+ x (0, T)) with 

IICIIl-(r+x(o,t)) < max{||Goid||L~(R+;iyi.p(n)), v^I^^I^}- (49) 
The relation (j48p now reads 

l|G||i=o(K+x(o,T);M'i.p(o)) < IICIIl°°(r+x(o,t)) exp(3Co||Vll||ii(o,T;Wi.p(n)))- (50) 

4.2.2 Estimate for the age derivate dgG 

We first remark that the derivative G' — dgG exactly satisfies the same PDE that G (see the first equation 
of ([57]) : that is due to the fact that u does not depend on the variable s). We then deduce the same kind of 
estimate that (|43)) : 

dt{\\G'\\LPin)) + ^a.(||G'|U.(o)) < CoHVM||^i..(o)||G'|Up(n). (51) 

But the initial conditions differ as follows: G'|(=o = 9sGoid and G'|s=o = SUeVlt. This last condition is 
obtained using s = in the equation (137^ . Note that this result is valid because we are working on regular 
solutions G„ (see the introduction of this proof) such that dtGn is continuous at s = 0. From the lemma [B. II 
given in the Appendix iBl we obtain for all (s, t) G R+ x (0, T) the estimate 

||G'(s,i,-)||LP(a) <C'(s,i)exp(Co^ II ViI|1m/i,p(o)) , (52) 
where C'(.,t)=<! " ^) if^<2ne., 
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For each s > 0, we estimate the L'^{0, T)-norni of the function 1 1-^ C'(s, t) as fohows: if T < Wc s then 



-2ne^'^lia.Goid||2p(o)(^')di' 



(53) 



< We ||i9sGoid||L^(R+;LP(n))- 



UT>Wes then 

rT /.aiTcs 



J C{s,trdt = J ||a.Goid|IL.(o)(s-^)di + 2Ue'-y \\Vu\\l,^^^^it - We s) dt 

= We ||9,Goid||L.(o)(i')di' + 2ne'- / \\Vu\\l,^^){t') dt' 

Jo Jo 



(54) 



< 21Je||a,Goid|ll.(K+;LP(o)) +2»e'' || Vm||2,.(o -^.^.(o))- 
Finally, we obtain (' e L°°(R+; L''(0, T)) with 

IIC'I1l~(M+;L'-(0,T)) < 2Ue-||asGold||L-(R+:LP(J7)) + 2Be|| Vm|1 L-(0,T;LP(f2)) • (55) 

The relation (j52p now reads 

l|G'||L«>(K+;L'-(o,r;LP(n))) < IIC'IU°°(R+;L'-(0,r)) exp(C*o||VM||Li(o,T;VKi'P(n)))- (56) 

4.2.3 Estimate for the time derivate dtG 

Isolating the term dtG in the equation ([37]) we have 



(57) 



< ^^I|G'||lp([2) + Co||M||vKi.p(n)ll'^G||iP(Q) + Co||G||wi,p(j2)||Vll||iP(f2). 



Introducing the Poincare inequality with a constant Cp = Cp(fl,p), holds since u vanishes on the boundary 
of the domain, we obtain 

WdtGh.in) < ^||G'|Up(o)+Co(l + Cp)||V«|Up(n)||G|l^i.i,p(n). (58) 
Taking the L''(0, T)-norm for the variable t, and next the L°°(M+)-norm for the variable s, we obtain 

|l(?tG||Loo(H+.i,.(o_T;LP(0))) <^llG'||i==(R+.ir(o,T;LP(a))) ^^g^ 

+ Co(l + Cp)||VM||i7-(o j-.iP(Q))||G||ioc.(R+x(o,T);iyi'P(0))- 

Using the previous estimates ([50)) and (|56)) . we deduce the result announced in the lemma l4?2l 

4.3 Estimates for the stress tensor r 
Lemma 4.3 Let T > 0, r e]l, +oo[ and p e]c?, +oo[. 

IfGe L°°(R+x (0,r); VFi'P(17)) anrf StG € L°°(R+; L''(0, T; LP(r2))) then the stress tensor r defined by 
the integral relation (1381) belongs to L°°(0,T]W^'P{n)) and its time derivative dtT belongs to L'^(0,T; LP(fl)). 
Moreover, there exists a continuous increasing function Fq : M+ i— > M+ such that 

||T||L~(o,T;Wi.p(n)) + \\dtr\\L'-{o,T;Lp{n)) < ^-Po(l|G||ioc(R+x(o,T);H'i'P(a)) + l|i9tG||Loo(R+.^r(o_T.iP(n))))- (60) 
The function Fq depends on fl, p and on the function growth of the function S. 
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Proof Since the function S is of class C^, we can introduce the following continuous and non-decreasing real 
functions 

5o : c e IR+ I — > max \S(G)\ G R+, 

\G\<c , ^ 

^ , ^ 61 

Si:ceR+^ max \S'(G)\ £ R+, 

|G|<c' 

where the derivative S'{G) denotes the 4-tensor whose the coefhcient {i,j,k,£), denoted di^ij)S{G)ki, is the 
derivative of (5(G)) ^.^ with respect to the tensor Eij of the canonical basis of the space £(]R'') of real d x c? 
matrices, see the Appendix \K\ 

Due to the continuous injection W^'^{Vl) ^ L°°(r2), the function G introduced in the hypothesis of the lemma 
is bounded in R+ x [0, T] x Q, and we have 

ll'5(G)||L'=o(R+x(0,T)xJi) < '5o(||G||ioo(R+x(o,T)xa)) 

(62) 

< 5o(Co |lG|j2;,c=o(R+x(0,T):lVi.P(n)))- 

To simplify, we note c := Co ||G|lioo(-B+x(o,T);Wi p(f2))- In the same way the function S'{G) is bounded by the 
real Si (c) . 

XP-norm for r — We easily have the following bound for the stress tensor t given by the formula ([55)) : 
[''■(^j x)\ < ^So{c) for a. e. {t, x) g (0, T)xV,. We note that we used m = 1. We deduce in particular that 

liT|lLoo(o,T;Lp(n)) < M''^So{c). (63) 
W^-'-'P-norm for t — Taking the spatial gradient of the expression ([55)) given the stress tensor we obtain 

VT(i, x) = — m(s)VG(s, i, a;) : S'{G{s, t, x)) ds. (64) 

2Ife Jo 

The meaning of the symbols here is the following. Component by component, the equality above written 

^ r+oo 

[VT(t,a;)]^^.^, ^diTjk{t,x) m{s)^d^Gi>m{s,t,x) d(^i„,)S{G{s,t,x)) ^^ds. 

poo 

Using the Holder inequality, the L^-bound on S'{G) and the fact that / m = 1, we obtain 

^0 

p POO 

|VT(i,a;)|P = — m{sf^VG{s,t,x):misy'^S'{G{s,t,x))di 
- ^'5i(c)P^°°m(s)|VG(s,t,a^)|f ds. 
Integrating for a; G f2, we obtain 

||Vr(i,.)||i,(,,) < ^SiicT J^^ m(s)\\WG{s,t,-)\\l,^^^ds. (67) 
Due to the definition of the bound c, this implies 

l|VT||ioo(o,T;Lp(n)) < ^'^i(^) (68) 

LP-norm for dtT — Similarly, we obtain a bound for dtT in LP{il): we have 

^+00 



(65) 



(66) 



/ m{s)dtG{s,t,x) ■.S'{G{s,t,x))ds. 



(69) 
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Using the Holder inequality to control the quantity \dtT{t, x)\p, and next an integration for x E ft, we obtain 
||5,r(t,.)||L(n) < ^S^icr J^^ m{s)\\dtG{s,t,-)\\l,^^^ds. (70) 
Due to the assumption on dtG which implies that c := H^f G||icxj(R+.ir(Q T^.^pfs^-))-) < +00, we deduce 

l|5t-r||i,-(o,T;Lp(a)) < ^'^i(^)^- (71) 

The estimates (1551) . and ((7T]) show that r and dtT are bounded respectively in L°°(0, T; W^'''(^^)) and in 
L^{0, T; LP{il.)), and that these bounds continuously depend on c and c, and increase in both variables. □ 

4.4 Proof of the theorem 13.11 

For any T > we introduce the Banach space 

^(T) = L'-(0, T; Wo'^in)) X C(K+ x [0, T];LP{n)) x C([0, T]; 2.^(17)) 
and for any i?i > 0, i?2 > and i?3 > the subset 

jr(r,i?i,i?2,i?3) - {{u,G,t) e ^(T) ; we L'-(0,r;i^(Ap)), dtueL'^{0,T;Hp), 
G e L°°(M+ X (0, T); W^1'P(17)), a^G, ^tG £ i°°(M+; i'^(0, T; X''(l}))), 
r e L°°(0, T; W^^'^ (f!)), a^r e ^''(0, T; (f])), 

M|t=o = uq, G\t=o = Gold, G\s=o = S, (72) 
ll'"llL'-(o,T;iy2-p(n)) + ||i9tM||L'-(o, T;LP(a)) < -Ri, 

||G|jLoo(u+x(0,T);H'i'P(0)) + |19sG|jL~(K+;L'-(0,T;LP(a))) + || C?tG|l l~ (K+ ;L'-(0,T;LP(n))) < -^2, 

Remark 4.1 S'mc/i a set is non empty, for instance if Ri and R2 are large enough. More precisely, if 

Ri > -. llwollDrCn) and R2 > ||Goid||L=(R+:Wi'P(0)) (73) 

1 — CJ p 

then for any T > and any R^ > we can build a velocity field u* such that (u* , Gold, 0) G J^(f{T, Ri, R2, R3), 
see an example of construction in fl5[ \21^ . 

Remark 4.2 // {u,G,t) e J^(T, _Ri, i?2, -R3) for some T, Ri, R2 and R3 then the velocity field u and the 
tensor G are continuous with respect to the time t and the age s. In fact, these continuity properties follow 
from the Sobolev injections of kind W^'°'{0, A; X) C C{[0, A]; X), holds for a> 1. Moreover, they make sense 
of the initial conditions u\t=o = Uq, G\t=a = Goid o-nd G\s=a = S. 

We consider the mapping 

$ : Jf{T,RuR2.Rz)~^ ^(T) 

{u,G,t) ^ {u,G,t), 
where u is the unique solution of the Stokes problem p6p with 

g = -meuVu+ dWT + f; (75) 

where G solves the problem (j37p and where t is given by the integral formula (|38| . The goal of this proof is 
to show that the application $ has a fixed point. For this we first prove that $ leaves a set J^{T, Ri, i?2, ^3) 
invariant (for a "good" choice of T, R2 and R3). 
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Let T > 0, i?i > 0, i?2 > 0, i?3 > and (m,G,t) G ^(T, i?2, i?3)- If we denote by (m,G,t) = $(M,G,r), 
we will show that the previous lemmas imply estimates of (it, with respect to the norms of (U, G,t), 

that is with respect to (T, i?2, ^s)- 

Velocity estimate — From the lemma HHI we can estimate u and dtU using the norm ||fl'||L' (o,T;LJ'(o)) ■ For 
the source term g given by the relation (1751) . we have 

\\9\\L^{0,T;LP{a)) < IIU • Vw||ir(o,T;LP(n)) + II''^IIl'-(0,T;W1'P(J7)) + II /II L'-(0,T;LP(a)) ■ (76) 
^ ^ ^ ^ ^ ^ 

Ti T2 

Since we have the bound R3 on r in L°°{0,T] M^i'P(fi)), the term T2 satisfies T2 < r-i?3. The bilinear term 
Ti is more difhcult to estimate. We follow the ideas of [15_ and we generalize their result to the d-dimensional 
case (the paper only deals with the case d = 3): 

Ti < ||M||L2'-(o,T;L~(n))||VM||i2T-(o_r;Lp(n)) 

_ _ (77) 

< ll«ll^2r.^^^^^^^(^^^||«IU-(0,T;Wi..(n)). 

But we have the following estimate (see pT]): 

p-d d 

INlL»(a) < C||w|U(j2)||m||^i,p(j^), (78) 
which, after integrating with respect to time, implies 

p~d d 

"'""L^(O.T;L°'(a)) - ^ll'"lli°"(0,T;LP(0))ll'"llL2.-(o,T;Wi.P(0))- C^^) 

Note that the constant C introduced here only depends on fl, p and d. Moreover, by interpolation, we have 

1 j_ 

l|w||L2.(o,T;Wi.p(n)) < ll'"lll~(o,T;Lp(n))ll"llL'-(o,T;vy2-p(n))- (80) 
Using ([75)) and (j80p . the estimate ([77)) now reads 

Ti < CT^|N|5!t,T;Lp(0))INl3fo,T;H'2,P(a))- (81) 
Finally, we use u{t, x) = Uo{x) + diu to obtain 

II^IIl°=(0,T;Lp(O)) < 111*01^^(12) + l|5tw||Li(o,T;LP(n)) 

< ll'"o||LP(f2) +r '■ ||(9tM||ir(o j-.ip(Q)). 

Using the bound Ri for It and diu given in the definition of the subspace ^(T, _R2, ^3), the estimate ([ST)) 
becomes 

r»-d P + '^ 3p-d 
R 2p 11^, J| 2p 

We now use this bound to control the source term g. The lemma HTT] implies 



Ti < CT^R^'" WuoWJ".^. + CT^-rRl (83) 



l|l'llL'-(0,T;H'2,p(n)) + l|c^tW||L'-(0,T;LP(n)) < ( II '"O || D'-(f2) + II /II L'-(0,T;LP(f2)) 

1 - w V (84) 
+ $neCr"i?f ||mo|II,j(i^) +^eCr*'i?2 +T^i?3), 

where the assumptions p > d and r > 1 imply that a, /3, 7 and 5 are positive numbers. This estimate 
can be rewrite as 



ll«IU-(o,T;iy2,P(o)) + ||attt|U,.(o,T;Lp(o)) < Ci (^1 + T" R3 + K{T, i?i)j , (85) 

where Ci may also depends on uj and on the norm of Uq and / in their spaces. It is important to notice that 
for each i?i > we have ]im^K{T, Ri) = 0, and that Ci > j^\\uo\\D^{n)- 
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Deformation gradient estimate — From the lemma we have 

l|G||ioo(R+x(o^T).VKi.p(f2)) + l|c^tG||Loo(R+.ir(n_T;Lp(n))) < ^2(1 + i?i)exp(C3T^"~i?i). (86) 
Stress tensor estimate — From the lemma we have 

\\'''\\L°°{0,T;W^-P{n)) + |19tx||ir(o,T;LP(n)) < 2^^o(^2)- (87) 

Invariant subset — If we then successively choose 
R* = 2Ci, 

i?2 = ^'2(1 + Rl)exp[C3Ri) + \\Goid\\ L°° (R+ -w^-p (n)) , 

and < 1 small enough to have t} R^ + K{T^,,Rl) <1, 



then we verify that J^(r*, i?*, i?2, ^^3) 7^ (that is the inequalities ([75)) hold). For such a choice, the esti- 
mates ([55]), dini) and dST]) imply that <I>(^(r-^, i?^, i?^)) C .yt^(T^, R^, Rt,, R3). Moreover the function $ is 
continuous and Jif{Ti,,Ri,R2,R^) is a convex compact subset of ^{Ti,), see [3T] for similar properties. We 
conclude the proof using the Schauder's theorem. □ 



5 Proof for the uniqueness result 



This section is devoted to the proof of the local existence theorem 13.21 

As usual, we take the difference of the two solutions indexed by 1 and 2. The vector u = Ui — U2, the scalar 
p = Pi — P2 and the tensors r = Ti — T2, G — Gi ~ G2 satisfy the following: 



' ^e{dtu + Ui ■ Vm + u ■ VM2) + Vp — (1 — a;) Alt = div -r, 
div u = 0, 

+00 

m(s) S{Gi{s,;-)) -S{G2{s,;-)) 



UJ 



ds, 



(89) 



1 



dtG + —dsG + Ml • VG + tt • VG2 = Gi • Vm + G • Vm2, 

together with zero initial conditions m|(_q — and G|^_q — G\^_^ — 0. Note that the regularity of Gi and the 
definition of the stress tensor implies that Ti £ L°°{0,T; W^''^{il)) (the proof is similar that those presented 
in the proof of the existence theorem 13. ip . The uniqueness proof consists in demonstrate that u = and that 
G = X — 0. We will initially provide estimates on these three quantities. 

Velocity estimate — Taking the scalar product of the first equation of the system ((89)) by u in _L^(r2), we 
obtain 



Die 



From the Holder inequality and the Young inequality, we obtain 

4 „ ,„ 



T ■ Vu — 5He / {u ■ VW2) • u. 



mcdt{\\u\\l2 



^n}) ^ (1 ^ ^)l|V'"llL2(a) < i_ ^ H' llL2(a) -'"^11 v"2||L°=(f2)||"|lL2(Q). 

Introducing Z{t) = *He ||tt||^2({2-) and Cz{t) = || Vm2||l=°(si) €E L^{0,T), this estimate reads 



-^4^U2\\L^(n)\\u\\l2 



Z'(i) + (l-a.)||VM||i.(,,) 



< 



I-LU' 



U^n)+Cz{t)Z{t). 



(90) 



(91) 



(92) 
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Stress tensor estimate — From the definition of the stress tensor r in the System ([89)) and using the 
Cauchy-Schwarz inequality, we have, for all {t,x) e (0,T) x fl: 



\Tit,x)\^< 



UJ 



(s) S{Gi{s,t,x)) - S{G2{s,t,x)) 



As. 



(93) 



By assumption, the function S is of class C^, so that S' is bounded on each compact. Since G^, i G {1,2}, 
belongs to L°°(R+ x (0,T) x i7) we deduce that there exists a constant C", only depending on the norm 
l|Gj||ioo(R+x(o,T)xn) such that |5(Gi) -5(G2)P < C'\Gi - G2P a.e. in M+x(0,T) x Vl. We deduce 

\r{t,x)\^<—^ m{s)\G{s,t,x)\^ds. 
yUe Jo 

Integrating with respect to cc G we obtain 



T|lL2(fl) < 



2ne' 



(94) 



(95) 



where we introduced Y{t) = I m{s) || G(s,t 

Jo 



2 ds. 



Deformation gradient estimate — Taking the scalar product of the last equation of the system ([5^ by G 
in L^(f2), we obtain 



^9*(||G||i.(^)) + ^9.(||G||i.(^)) 



/ (Gi • Vm) • G + / (G • VM2) G [ {u- VG2) • G. (96) 

Jf2 Jo JO 



Using the Holder inequality, we have the estimate 



1 



1 



2^^t(l|G|li.(o)) + ^5.(||G||i.(,,)) <||Gi|U»(o)||V«|U2(s,)||G|U. 



(O) 



iVtt2||L~(0)||G||i. 



(O) 



(97) 



lM|l^^^^^||VG2|U-(a)||G|U2(n). 



Due to the Sobolev continuous injection H^{il) ^ L<i-2 ($7), the Poincare inequality and the Young inequality, 
we obtain for all e > 0: 



1 



+ ^a.(||G||i.(^)) < e||Vix||i.(j,) + CrmGWh 



(98) 



where the function Cyit) - snp,^^^ {f l!G'i||i=o(n) + I1Vm2|1l-(0) + ^ I! VG2|l|.(j^) } G L^{0,T), and where 
the constant C depends on fl, p and d. Multiplying this estimate (j98p by to(s) and integrating for s G (0, +00) 
we obtain 



Y'it) + ^J^ ™(5)5.(||G||i.(^))ds<e||VM||i.(f,)+Cy(i)r(i). 



(99) 



Using a integration by part, the integral becomes 

f +00 







-m'(5)||G||i.(o)ds+ m(s)|lG||i.(f,) 



Using the following arguments: 

- The memory function m is non- increasing, that is — m' > 0; 

- The function G is bounded on M+ x (0,T) xfl and m G L^(M+), that is lim m{s)\\G\\'i2,m 

- We have the following development with respect to the variable s for G: 



= 0; 



(100) 



G(s) = G| ^ + sdsG\ „+o(s) sWeVu. 

^ ' I s— Is— ^ ' 
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Moreover, m € L^(R+) so that liin to(s)|1G|1|,2(q) = 0. 

Hence the integral is non-negative so that the estimate ([M)) now reads 

Y'{t)<e\\Vu\\l,^,,)+CY{t)Y{t). (101) 

Uniqueness result — Finally, adding and (|10ip with the choice e = 1 — w, and using the estimate (IMI) . 
we obtain 

(r + z)'(t) <Cyz(t)(r + z)(t), (102) 

where the function C-yz is a linear combination of Cy and Cz- In particular we have Cyz € L"'^(0,T). The 
Gronwall lemma and the initial condition 1^(0) = Z{0) = imply that Y = Z = 0. We deduce that m = and 
G = and that consequently the stress r = and the pressure p is constant in (0, T) x il. □ 

6 Proof for the global existence with small data 

This section is devoted to the proof of the local existence theorem 13.31 Arguing as in the proof of the 
theorem l3. 11 we introduce the space ^(T), the subspaces J^{T, Ri, R2, R3) and the mapping $. 
For {u, G,t) S J^{T, i?2, R3) and (it, G, x) = ^{u, G,T) recall that we have the following estimates (see 
the estimates (EH), (EH) and ((57)) ): 

C'l ( 

ll'"IU'-(0,T;W2.p(O)) + ||c'tw||L'-(0,T;LP(O)) < ( II '"O II + I1/I1l'-(0.T;Lp(O)) 

1-wV (-^Qg^) 

+ $HeCr"i?f ||Mo||2,.(i^) + IRe Cr'^i?2 + T^i?3^ 

l|G'||ioo(K+x(o,T);H'i.p(n)) + llf^iG'||Loo(R+.ir(o,T;Lp(a))) < (^2(1 + -Ri)exp(C3r^^"i?i). (104) 

II''"||l°°(o,T;H'1.i'(o)) + l|5t-r||ir(o,T;Lp(n)) < 25^'^°''^2)- (105) 

Note that the constants Ci , C2 , C3 , C and the function Fq introduced in these three estimates do not depend 
on uj. For a time T > given, we successively choose 

/?* - ~ ^ 
^ ^ 2Ci 9^0 6^-5' 

i?^ = C2(l + i?t)exp(C3Ti-^i?|) + ||Goid||L~(R+;Wi-(n)), (l'^^) 
2Be" 



^3 - ^J^^RD- 



We verify that, for Mqi / small enough in their norms, and for w small enough too, J^^(T, i?*, i?2i ^3) 7^ ^ (that 
is ([75]) holds). For such choices oi R*, R2, to and small norms of Mq, / and Gold, the inequalities (I103p . (jl04p 
and (fT05)) imply that <I>(Jf'(r, i?!, i?^, i?^)) C .^{T,R\,R*2,R\). Moreover the function $ is continuous and 
,ytf{T, R*, i?2, ^3) is a convex compact subset of 3§{T). We conclude the proof using the Schauder's theorem 
again. □ 

7 Proof on the stationary problem 

The stationary problem associated to the system (IT51) reads 

' <ne M • Vm + Vp - (1 - uj)Au = div T + /, 
divM = 0, 

r{x) ^ ^ m{s) 5(G(., x)) ds, (107) 

r^dsG + uVG^GVu, 
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coupled with the "boundary" conditions u\g^ = and = S. 

We win introduce a Banach space a convex compact subset J(f{Ri,R^) and a continuous mapping $ : 
,J^{Ri,Rz) — > in such a way that the system (|107l) is equivalent to a fixed point equation for $ in 
jr(i?i,i?3). More precisely, let us set = x LP(f7) and 

Il'"||w/2,p(f2) < ||T||vi/i.P(n) < ^3 |- 

For each i?i > and i?3 > it is obvious that we have J{f{Ri, R3) 7^ 0. We now consider the application 
$ : {u,t) £ J^{Ri, R^) I — > {u,t) E where u is the unique solution of the stationary Stokes problem 

divM = 0, (109) 

with the following source term g = — fHelt • VU + div r + /; and where t is given by the integral formula 

r{x) = ^ ™(s) 5(G(s, x)) ds, (110) 

the tensor G being the unique solution of the equation 

-^dsG + u-WG = G-Wu, 

Wc (111) 
G\ ^^5. 

Remark 7.1 In the proof of the theorem \3.4\ for reasons of clarity, we often only indicate the dependence of 
the constants with respect to the parameter u). We write A < B if there exists a constant C which does not 
depend on uj (but depending possibly on fJ, p, d,...) such that A < CB. 

7.1 Estimates for the velocity field u 

It is well known (the proof is given in [28]) that the unique solution of the Stokes problem p09p belongs in 
W^ P{n) as soon as the source term g G LP{U). Moreover we have the following continuity estimate: 

||VM||vyi,p(o) < -7-^\\g\\LPin), (112) 

where the constant C4 only depends on the domain and the integer p. For our purposes, the expression of 
the source term g allows us to deduce 

l|VM||vvi,P(n) < Y^(ll/IUno)+^e||VM|rvFi,P(n) + ll^llH'i-(o))- (113) 

7.2 Estimates for the deformation gradient G 

Following the same method that those to obtain an estimate of G for the unsteady case (see the subsec- 
tion we prove that the solution G of the system (|llip (which exists and is unique by the characteristic 
method) satisfies 

— 4(||G||wi-p(o)) < 3Co||Vu||M/i,p(o)|lG|jH'i.p(n), (114) 

where Co = Co{fl,p) corresponds to a constant of the continuous injection W^'P{il.) ^ L°°(ri). Since 
IIG'llvKi'P(a) |s=o = V^l^^l^, the classical Gronwall lemma implies that for all s > 0, 

l|G(s,-)||wi-(n) < Vd\n\hxp(3CoWe\\Vu\\w^.,(^n}s). (115) 
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To give sense to the "initial" condition G\s=o — S we need estimate on dsG. Isolating this term on the 
equation (jllip we obtain 

\\dsGis,-)\\L.^n)<'rBc[\\Gis,-)\\w^.,^n)\\'^^\Lnn)+ (116) 
That is, due to the estimate (|115|) . for all s e M+: 

||a,G(s,.)||LP(o) < ||VlI||iP(o)exp(3Co2ne||VM||M/i.P(n)s). (117) 
7.3 Estimates for the stress tensor r 

To control the stress tensor given by the integral (jllOp . we must control the quantity S{G). Using the 
assumptions (A4) on the growth of this function: \S{G)\ < \G\°-, \S'{G)\ < \G\'' we deduce that for all x e Q: 

/"OO 

\t{x)\p<ujP m{s)\G{s,x)rds, (118) 
Jo 

(119) 



and 



POO 

\S/t{x)\p<ujP / m{s)\VG{s,x)\P\G{s,x)\''Pds. 
Jo 

Since G G W^'P{i}) C L°°(f2) (with respect to the spatial variable), we have 

/•OO 

MUn) ^ m{s) \\G{s, Oll^^.^c^) d., (120) 

and 

poo 

W^^WUnj^^'l m{s)\\G{s,-)C\'JP^^ds. (121) 
Using the estimate (|115p on G, and the assumption (A3), that is 7n{s) < e~"^, we obtain 

\\T\\w^'P(n) < exp((3Co2IJccp - q;)s) ds, (122) 

Jo 

where c = max{a, 6+1} and where C5 does not depend on uj. 

7.4 Proof of the existence result for the theorem 13.41 

Let (u,t) e J{f{Ri,Rz) and {u,t) = $(It,T). Using the relations (|113p and (|122p . we have 

I|Vm||wi,p(o) < Y^(ll/IUno)+^e^? + ^3), (123) 

pOO 

Jo 

Taking i?i smah enough to have 3Co2Uccpi?i < a/2, the inequality (|124p becomes 

lk|k.,.(o) < (125) 

We now choose — ^"'^■'^ and we verify that for / and to small enough, the relations p23p and (jl25p 
imply the inclusion ^{,J€{Ri,Rj)) C ,J^{Ri,R-i). As for the proof of the theorem 13. 11 the Schauder'theorem 



implies the existence of a solution with estimates for u and for t. Next, using (|115p and (jll7p . we deduce the 
estimates for G and for dsG. 
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7.5 Proof of the uniqueness result for the theorem 13.41 

Proceeding as in the proof of the theorem 13.21 (see section[5]), we consider two solutions {ui,Ti) £ J^{Ri,R3) 
with i S {1,2}. We note Gi the corresponding deformation gradients. From the estimate (IllSp . we have for 
all s e M+: 

\\G,{s,-)\\wi.,^n) < Vd\n\h'^'<'^''^^'. (126) 

We next introduce the differences u = Ui — U2, t = ti — T2 and G = Gi — G2 and consider the equations 
satisfied by these differences. 

• The scalar product of the Stokes equation for tt by m in L^{il) gives 

(l-^)l|Vw||i.(n) < ^^||r||i.(^)+^Re||Vtt2||L»(o)||M|li2(j,). (127) 

• From the definition of x as a difference of two integrals depending on S{Gi) and S{G2) we obtain 

\t{x)\<—J^ mis)\f'{G,{s,x))\\Gis,x)\ds, (128) 
where WG^is, ■)\\w^..in) < max {\\G,{s, ■)\\w^,,^n)} < e3C«2"^«l^ Using the bound |/'(G3)| < {G^l' we 

iG{l,2} 

obtain 

/•OO 

\t{x)\<uj m{s)\G3{s,x)\'\G{s,x)\ds. (129) 

JQ 

By the continuous injection W^'P{fl) ^ L°°{n), and from the estimate on ||Gi||iyi.p(o)i * G {1,2, 3}, we have 

poo 

\t{x)\<uj m{s)e^^°^'^''''\G{s,x)\ds. (130) 
Taking the square and integrating with respect to x (after using the Holder inequality) we deduce 

Mhm<u^' / m(.)e6'^°™»«^''^||G(s,.)|li.(o)ds- (131) 
Jo 

• Taking the scalar product of the equation on G by G in L^(fi), we obtain (see (|n7| for similar calculations) 
-4(||G||i.(jj)) <||Gi|U^(jj)||VM|U2(f,)||G|U2 
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+ l|VM2||L-(f2)||G|||2(t2) + ||VG2||Ld(o)||VM||i2(o)||G||t2(n). 

From the Young inequality, this relation is written 

4(||G||i.(j,)) <A(s)||VM|li.(^)+i3|lG|!i.(^), (133) 

where A{s) = ^(HGiH^^^jj^ + II VG2|l^d(f2)) and B = 2 2ne(|| VM2||Loo(n) + -Ri). By the Gronwah lemma 
(and using the fact that G{0,x) — 0) we obtain 

||G(s, •)||i2(n) < /' e^^^'A(s')|| Vwlli.(^) ds'. (134) 
Jo 

Using the bound on G, we have A{s') < .^e^'^o^'^^i* for ah s' e M+. We deduce 

l|G(.s,-)lli2(,,) < ^(^%(6^o=«eii.-BKd,.)|,^„||2^^^^^ (135) 
Now, the estimate (jl31l) becomes 

^ ( m(s) e(6Co2Hcflifc+B). g(6Coa»ci?i-B).'d^/ ^ II VM||i.(f,). (136) 



2 < 
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Assuming (A3), that is m{s) < e the integral satisfies 

Hence, for Ri small enough (such that 2Wz{3Cob + 2)Ri < a, we will note that B < 42Uei?i) the integration 
with respect to the variable s in the last integral converges and we have 

1 r°° 

a< . / p(6Co21Tcfli(fc+l)-a)s' - / /-, oo^ 

^ a-(6Co2IIei?i6 + i3)yo ^ ^ 

Finally, if in addition we choose i?i small enough to have 6Co2ITei?i (^^ + 1) < a then the integral J" converges 
and we have the following bound for the stress: 

lklli^(o)<^l|V«|li^(n)- (139) 

Using the relations (jl27p and (I139p . we deduce that for R\ and lo small enough we have m = and r = 0. 
Next, the relation (I134p implies G = 0, that concludes the uniqueness proof of the theorem 13.41 □ 



8 Conclusion 

In this article we are interested in the mathematical properties of models of viscoelastic flow. We have shown 
that many known results for differential laws could be adapted to integral models. Nevertheless some differences 
persist: 

• Our results are formulated in the U —LP context (following [H]). It seems possible to reformulate them 
for more regular solutions in a H" context (following [21]). 

• Some results: the global existence with small data (theorem l3.3p and the stationary study (theorem l3.4p . 
are proved for the relaxation parameter cj small enough only, that is to say for the flows which are not too 
elastic. About the differential models, this assumption can be removed, see for instance [TUl 130] for the 
global existence, and 15 for the stationary problem. In these two cases, the results about differential 
models strongly use the structure of the equation, and it seems difhcult to adapt such methods for 
integral models (see also the remark [3TT]) . 

• There exists differential models which have no apparent equivalent in terms of integral models, for 
instance the co-rotational Oldroyd model. This study does not cover these cases. Similarly, there are 
also integral models more general than those studied here. In 36 , R. I. Tanner introduce models where 
the memory m also depends on the invariants 1\ and I2 ■ It might be interesting to study from a theoretical 
point of view these models and to observe if the approach taken here fits well. 

On the other hand, it is possible that classical integral models perform better than differential models of 
Oldroyd type. In fact, most of these models have a stress which is naturally bounded via the definition of S, 
see the examples given by the equation (PU)) . (PT]) or ([5^ . and the Appendix \X\ While getting weak solution 
seems very difficult, knowing a priori bound on the stress is an interesting information (see [9 for an example 
of a criteria for the explosion in the Oldroyd model) . 

Finally, the theoretical results shown in this paper allow to consider a lot of work on models of integral type. 
The well-known results for the differential model can be generalized to integral models. For instance, the one 
dimensional shearing motions and Poiseuille flows admit global existence for usual differential models, see [19]. 
In this regards, the works of A.C.T. Aarts and A.A.F. van de Ven [T] are interesting: they study the Poiseuille 
flow of a K-BKZ model. Would it be possible to prove global existence for such one dimensional flows when 
we use more general integral models ? An other possible generalization concerns the behavior of viscoelastic 
flows in thin geometries (in the fields such as polymer extrusion or lubrication), or in thin free-surface flows 
(for study mudslide or oil slick). Recent works [Sj [2] and |7] can provide answers to the differential models, 
and we can imagine the same kind of works for integral models. 
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A Tensors and the strain measure function S 

A.l Some remarks on the invariants of the Finger tensor (d = 3) 

For a matrix B e £(R'^), we usually define three invariants: 

/i=Tr(i?), /2 = i((Tr(B))2-Tr(B2)), /g^detS. 

We specify in this subsection some properties of these invariants in the context studied here, that is when B 
represents a Finger tensor of an incompressible flow. 

First of all, this incompressibility condition implies that detS = 1. Consequently the third invariant I3 is 
useless. Next, using the Cayley-Hamilton theorem we have B^ = B^ — IiB + I2S and we deduce that 

I2 = Tr(B-i). 

By definition, B = '^F ■ F is real positive-definite matrix, and consequently it is diagonalizable. Using a basis 
formed by its eigenvectors, we have Ii — Xi + \2 + A3 whereas A1A2A3 = 1. An inequality of arithmetic and 
geometric means indicates that /i > 3, and in a similar way we prove that /2 > 3. We deduce that for all 
/3 e [0, 1] we have 

(3h + (1 - /3)/2 > 3. 
That mathematically justifies the PSM model ([50]) and the Wagner model (PT|) . 

A. 2 Notion of derivative for the Strain measure tensor 

The strain measure function S : £(R'^) — > C(R'^) can be viewed as an application M.'^^'^ — > R''^''. For 
^ ^ i, j ^ d, we introduce the matrix 

Eij = 



as element of the basis of C(W^) and we define the differential of S from its Jacobian S' :~ f-^^r—l . Note 

dSki 

that in this paper, we use the notation d^ij-^Ski := . 

oEij 

A. 3 Norm for the 4-tensor 

The notion of derivative introduced above involves the use of tensors of order 4. Recall that for a 2-tensor 
A = {A)ij, we use the usual algebra norm defined by |Gp := Tr(-^G • G). For a 4-tensors H ~ {H)ijki, we 
introduce the following algebra norm: 

We will note that this norm having the following property: |A ® B| = \A\ \B\ for any 2-tensors A and B. 
A. 4 Example for a PSM model 

Consider the example function corresponding to a PSM model (see the equation ([50)1 with a — A and (3 ~ 1): 

S:Ge C(R'^) ^ ^— — e C(R'^) where B=^G G. (140) 

^ ' 1 Tr(B) ^ ' ^ ' 




Proposition A.l For all G e C{R'^) we have \S{G)\ < 1. 
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Proof Using the norm on the 2-tensor, we have for all G E C{W^): 

that implies that S is bounded by the constant 1. □ 
Proposition A.2 For all G e we have \S'{G)\ < 

Proof The derivative of the function 5 is a function with values in the set of 4-tensors: 

^ %)(Bfc^) %)(Tr(B))Bfc; 

l + Tr(B) (l + Tr(B))2 
_ {^kjGii + SijGik) 2GijBki 



(142) 



l + Tr(B) (l + Tr(B))2' 

Taking the 4-tensor norm, we deduce that 

2|,5||G| 2|G|3 

''^ ^'^^'-TTW (1 + |GP)2 (^^^^ 

Since |5| = V^, this implies that \G\ \S'{G)\ is bounded by 2(1 + V^). □ 

B Gronwall type lemma 

Lemma B.l Let f : i— > a positive and locally integrable function. If a function y : x M+ i— > M 
satisfies, for all (s, i) e M+ x R+; 

5ty(s, + ^a,2/(s, t) < fit) y{s, t) (144) 
then we have, for all (s,i) g R+ x ]R+; 

y{s, t) < C(s, exp ^ /(t') di'^ , (145) 

. ^ iyis-:^,0) ift<Wcs, 
where C{s,t) = < ^ We' ^ ^ ' 

[ 2/(0, sues) ift>Wes. 

Proof Introducing the new variables u — ^ (SUe s + t) and v = ^ (2ITe s — i) , we can write the first equation of 
the system (|144p as a system on the function z{u, v) = y{s, t): 

duz{u,v) < f{u — v) z(u,v). (146) 

Since the function / is locally integrable, we obtain 



du 



z{u,v)eM - / fit')dt' 







< 0. (147) 



Integrating this relation between \v\ and u, we deduce 

z{u,v)exp(^~ f{t')dt'^ <z{\vlv)e^p(^~ J^^ /(t')dt')- (148) 

Due to the positivity of the function /, the exponential term in the last equation being less than 1. According 
to the sign of v, we have z(|w|,u) = y{0,t — QUe s) or z(|i;|,t)) = y{s — 357,0). That implies the result (|145p 
announced in the lemma. □ 
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